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20  ABSTRACT  (Cm Mw*2) 

Many  orographic  affects  are  produced  similar  to  that  obeerved  for  typhoons  passing 
over  Taiwan.  The  storm  tends  to  translate  at  about  twice  the  speed  of  the  basic  flow 
near  the  mountain,  while  its  intensity  is  reduced.  Air  flows  mostly  around  the  mountain 
range  instead  of  over  it,  forming  a  ridge  on  the  windside  and  a  trough  on  the  leeside 
slopes.  The  disturbed  flow  produced  a  mean  cyclonic  circulation  around  the  mountain. 
Therefore,  the  model  tropical  cyclone  makes  a  cyclonic  curvature  in  its  path  around  the 
north  end  of  the  island  mountain  in  its  path. 

Further  numerical  experiments  demonstrate  that  cumulus  heating  forces  the  cyclonic 
circulation  around  the  mountain.  We  found  that  the  lower  level  circulation  of  an 
unforced,  quasibaro tropic  vortex  is  blocked  by  the  mountain  range.  Instead,  secondary 
low-level  circulation  center  forms  in  the  induced  lee  trough.  The  secondary  low-level 
center  develops  as  die  upper  level  center  becomes  a  phase. 

A  vorticity  budget  is  performed  for  the  700-mb  airflow  just  before  the  landfall.  It 
confirms  the  importance  of  diabetic  processes  in  producing  the  obeerved  orographic 
effects.  Diabetic  processes  generate  convergence  to  maintain  the  vorticity  of  the 
tropical  cyclone.  The  horizontal  advection  of  positive  vorticity  in  conjunction  with  the 
leeside  vortex  stretching,  results  in  the  mean  positive  vorticity  around  the  topography. 
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THE  OROGRAPHIC  EFFECTS  INDUCED  BY  AN  ISLAND  MOUNTAIN  RANGE 
ON  PROPAGATING  TROPICAL  CYCLONES 


1.  INTRODUCTION 

The  island  of  Taiwan,  with  its  Central  Mountain  Range  (CMR) 
soaring  up  to  more  than  3000  m  at  only  50  km  from  its  eastern  sea 
shore  (Fig.  1),  Induces  many  interesting  topographic  effects  on 
the  weather.  In  winter,  the  northeasterly  monsoon  winds  travel 
unobstructedly  over  the  kuroshlo  and  produce  prolonged  periods  of 
rain  over  most  of  the  northern  and  eastern  parts  of  the  island. 

In  summer,  the  northwestern  Pacific  high  pressure  system  is  in 
control  and  the  weather  is  highlighted  by  frequent  typhoons.  As 
typhoons  approach  the  island,  their  behavior  can  be  influenced 
substantially  by  the  high  CMR. 

Brand  and  Blelloch  (1974)  found  that  westward  traveling 
typhoons  tended  to  move  cyclonically  around  the  northern  side  of  the 
CMR  (Fig.  1).  The  average  translation  speed  changed  from  8  m  s"*  48  h 
before  the  landfall  to  12  m  s'**  at  the  landfall.  Noticeable  weak¬ 
ening  started  12  h  prior  to  the  landfall.  Chu,  Wang  and  Pao  (1977)  showed 
by  streamline  analysis  some  detailed  flow  structures  such  as  flow 
separation  and  vortex  shedding  associated  with  the  typhoon  circula¬ 
tion  around  the  island.  Wang  (1980)  has  compiled  so  far  the  most 
complete  data  of  53  typhoons  that  came  close  to  Taiwan  between 

1946-1975.  He  documented  of  these  typhoons  their  paths,  intensities, 
propagating  speeds  and  life  histories.  He  found  that  some  typhoons 

passed  over  the  island  by  taking  a  cyclonic  track  similar  to  those 

described  in  Brand  and  Blelloch.  There  were,  however,  some  other 
cases  where  typhoons  passed  over  the  island  by  forming  a  secondary 

circulation  center  to  the  leeslde  so  that  the  storm  tracks  were  not 
continuous. 

Muraacript  submitted  S*ptemb«T  16. 1981.  . 


Th«  observed  behavior  of  typhoons  has  been  reproduced  in  laboratory 
experiments  (Chang  and  Chen,  1969;  Pao  and  Hwang,  1977;  Wal,  1979). 
Unfortunately,  neither  quantitative  analyses  nor  physical  insights 
in  understanding  the  mechanism  of  the  island-induced  effects  were 
given. 

This  study  is  to  Investigate  the  physical  mechanisms  Involved 
in  the  observed  behavior  by  numerical  simulations  with  an  idealised 
topography.  We  will  discuss  in  the  following  Sections  2  and  3, 
respectively,  the  numerical  model  and  the  incorporation  of  the 
idealized  topography  to  the  model.  The  orographic  effects  of  the 
island  mountain  range  on  an  approaching  tropical  cyclone  will 
be  discussed  in  Section  4,  We  will  examine  in  Sections  5  and  6,  the 
effects  of  cumulus  convection  and  surface  friction  on  tropical  cyclones 
moving  over  a  mountain  range.  A  vorticity  budget  will  be  given  in 
Section  7  to  explain  the  observed  cyclonic  movement  of  the  typhoons 
near  the  CMR. 

2.  THE  NUMERICAL  MODEL 

The  problem  of  the  air  flow  over  an  Isolated  mountain  is  a 
very  Interesting  one  in  fluid  dynamics  as  well  as  in  weather  fore¬ 
casting.  In  earlier  mathematical  models  (see  Smith,  1979  for  review), 
the  case  of  stratified,  rotating,  barotropic  basic  fluid  over 
mountain,  whose  surface  coincides  with  an  isentroplc  surface,  is 
considered.  The  solution  of  such  a  quasi-geos trophic  approach  is 
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that  the  flow  near  the  mountain  la  similar  In  Its  structure  to  the 
classical  description  of  a  cold-core  anticyclone.  The  maximum 
vertical  displacement  of  the  lsentroplc  surfaces  decreases  with 
height,  but  the  horizontal  extent  of  the  lifting  increases  with 
height  (  see  e.  g.  Bussl  and  Tibaldi,  1977).  Thus,  quite  different 
from  the  simple  conservation  of  potential  vortlcity,  an  air  column 
approaching  the  mountain  would  undergo  a  stretch  before  the  moun¬ 
tain  elevation  due  to  the  lifting  of  the  stream  surfaces,  a 
contraction  over  the  mountain,  and  then  a  stretch  after  passing 
the  mountain. 

It  is  very  difficult  to  find  a  meteorological  situation  that 
fits  the  foregoing  discussion,  because  those  solutions  are  pertinent 
only  to  large  scale  (~  1000  km)  mountains.  For  the  mesoscale 
mountains  under  consideration  where  the  Rossby  number  (Ro)  approaches 
1,  the  quasi-geo8trophic  assumption  breaks  down  and  inertial  effects 
must  be  included.  Merkine  and  Kalnay-Rivas  (1976)  used  the  semi- 
geostrophic  approximation  (Merkine,  1975)  to  study  the  rotsting 
stratified  flow  around  finite,  isolated  topography.  They  found  the 
flow  consists  of  a  topographically-bounded  anticydonlc  vortex.  The 
vortex  is  baroclinlc  near  the  mountain  and  barotroplc  far  away  from 
the  mountain.  The  anticydonlc  vortex  intensifies  the  velocity 
field  on  the  left-hand  side  of  the  mountain  (facing  downwind)  and 
decreases  it  on  the  right  hand  side.  They  also  found  from 
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trajectory  analysis  that  all  trajectories  possess  cyclonic  curvature 
over  the  mountain.  The  cyclonic  curvature  must  be  accompanied  at 
higher  elevation  by  anticyclonic  shear  in  order  to  be  consistent  with 
the  conservation  of  potential  vorticity.  Buzzi  and  Tibaldi  (1977) 
obtained  a  similar  solution  using  an  expansion  in  powers  of  Ro  to 
solve  for  the  flow  in  an  unbounded  fluid.  The  semi-geostrophic 
solution  remains  symmetric  about  the  topographic  axis  perpendicular 
to  the  flow.  Bussi  and  Tibaldi  stated  that  the  semi-geostrophic 
assumption  is  still  too  restrictive  to  describe  tne  propagation 
of  inertial  gravity  waves,  especially  when  Ro  increases  to  or 
exceeds  unity. 

Hayes  and  Williams  (1977)  used  an  inviscid,  adiabatic,  primi¬ 
tive  equation  (PE)  model  to  study  the  steady,  unsheared  flow  past 
a  finite  large  scale  mountain  range.  Their  results  indicated  that 
the  mountain  induces  a  pressure  ridge  over  the  mountain  and  a  low 
level  trough  on  the  leeside.  They  showed  further  that  the  results 
are  insensitive  to  the  vertical  resolution  of  the  model  if  the 
horizontal  resolution  is  sufficient.  The  horizontal  resolution 
affects  the  fine  scale  features  of  the  flow.  Recently,  Tibaldi, 

Bussi  and  Malguzzi  (1980)  investigated  the  orographic  effects  of 
the  Alps  on  an  approaching  barocllnic  wave  with  a  PE  model. 

They  found  that  a  small  scale  barocllnic  process  caused  an  ampli¬ 
fication  of  the  initial  disturbance  produced  by  the  mountain  when 
the  large  scale  wave  passed  over  it. 
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In  our  present  study  of  a  diabatieally-forced  tropical 
cyclone  interacting  with  a  mesoscale  mountain,  an  initial  value, 

PE  model  seems  most  appropriate.  The  numerical  model  we  used  is 
identical  to  the  one  in  Chang  and  Madala  (1980). 

a*  Governing  Equations 

The  governing  equations 

for  conservations  of  momentum, 

normalized  pressure  p/p  (=  a) 

s 

system  is  assumed  hydrostatic, 

b.  Physical  Parameterizations 

A  boundary  layer  parameterization  based  on  a  generalized 
similarity  theory  (Chang,  1981)  is  employed.  The  surface  layer 
logarithmic-linear  profiles  are  matched  with  the  mixed  layer  profiles 
so  that  only  one  model  layer  is  needed  to  represent  the  boundary 
layer.  The  relevant  universal  functions  of  stability  A,  B,  C  and  D 
are  formulated  according  to  Yamada  (1976).  The  surface  roughness 
length  over  the  ocean  surface  are  computed  following  Charnock's 
equation, 

2  -  c  /  g  (1) 

o 


are  in  pressure-weighted  flux  form 
temperature,  and  water  vapor.  The 
is  the  vertical  coordinate.  The 
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The  roughness  over  the  mountain  is  assumed  to  be  a  linear  function 
of  elevation 

ZQ  m  1  +  0.0005  x  hg  cm  (2) 

We  note  that  the  PBL  over  rugged  mountain  surface  under  high 
wind  conditions  may  not  be  well  modeled  by  our  parameterization.  One 
shortcoming  is  the  assumption  of  constant  PBL  height.  The  PBL 
parameterication  is  nevertheless  applied  in  the  expectation  that  the 
similarity  theory  in  general  will  still  hold  under  these  conditions. 

A  simplified  version  of  Anthes'  (1977)  cumulus  parameterization 
is  utilized.  The  total  latent  heating  is  determined  by  vertical 
integration  of  water  vapor  convergence,  which  occurs  mainly  in  the 
boundary  layer.  Vertical  distribution  of  latent  heating  depends  on 
the  conditional  instability.  The  vertically  averaged  relative  humidity 
(RH)  regulates  the  partition  of  latent  heating  and  moistening.  Large 
scale  ascending  motion  in  a  saturated  environment  results  in  the 
release  of  latent  heat  on  the  resolvable  scale. 

c.  Initial  Conditions  and  Boundary  Conditions 

The  mean  hurricane  sounding  (Sheets,  1969)  is  used  for  the 

initial  thermodynamic  state.  The  lapse  rate  is  conditionally 


unstable  from  the  surface  to  about  350  mb  and  the  RH  Is  high  up  to 
500  mb.  The  Initial  flow  field  Is  a  uniform,  unsheared  easterly 
current  at  a  speed  of  5  m  s“^. 

The  initial  vortex  embedded  in  the  easterly  flow  is  non- 
divergent  with  tangential  velocities  described  by 

vo  (r)  “  vmax  ip  exp.  {  [  1  -  (-=-  )2  ]  }  (3) 

lmax  max 

where  the  maximum  tangential  velocity  vmay  is  16  m  s-^  at  a  radius 
of  rmax  *  210  km  from  the  storm  center.  The  surface  pressure  and 
temperature  fields  are  obtained  by  a  non-linear,  non-divergent 
initialization  method,  and  are  in  gradient  balance  with  the  flow 
field. 

The  Neumann  lateral  boundary  conditions  are  applied  to  all 
variables.  Sponge  boundary  conditions  are  imposed  at  the  inflow  and 

outflow  boundaries  to  prevent  bounding  reflection.  The  non- 

-2  -3 

dimensional  diffusion  coefficients  :<H  (At)  (Ax)  there  are  -''0(l6~-3) , 
which  are  one  order  of  magnitude  larger  than  those  in  the  interior. 
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d.  Model  Structure 

The  model  atmosphere  from  Pg  to  P  ■  0  is  divided  into  seven 
sigma  layers  as  in  Chang  and  Madala  (1980).  The  sigma  layers  are 
bounded  by  a  ■  0,  0.1,  0.2,  0.35,  0.55,  0.77,  0.93,  and  1.  All 

variables  except  a  ,the  vertical  velocity, are  defined  at.  the  center 
of  each  layer  and  0  is  defined  at  the  boundary  of  each  layer. 

The  momentum  points  and  mass  points  are  fully  staggered  in  horizontal 
direction  following  the  scheme  C  in  Arakawa  and  Lamb  (1981). 

The  horizontal  resolution  is  uniformly  60  km  throughout  the 
model  domain.  The  model  domain  is  3000  x  3000  km  (51  x  51  grid 
points)  so  that  the  movement  of  the  vortex  is  not  affected  by  the 
lateral  boundaries  over  the  period  of  integration. 

The  spatial  differencing  is  second-order  accurate  and  conserves 
mass,  mementum,  and  enthalpy.  The  split-explicit  method  (Madala,  1981 ) 
is  applied  for  temporal  integration. 

3.  The  Idealized  Topography 

The  idealized  topography  (Fig.  2)  is  used  to  represent  the 
terrain  of  Taiwan.  Due  to  the  horizontal  resolution  of  the  model,  some 
of  the  details  in  the  topography  shown  in  Fig.  1  are  smoothed  out.  The 
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idealized  topography  has  an  east-vest  extent  of  240  km  and 
north-south  extent  of  480  km  (Fig.  2a).  Fig.  2  also  shows  the 
east-west  cross-section  AB  (Fig.  2b)  and  north-south  cross  section 
CD  (Fig.  2c).  The  idealized  terrain  is  symmetric  about  its  ridge. 

The  topography  is  introduced  by  a  linear  growth  from  0-12  h, 
because  a  dynamically  balanced  initial  state  with  the  orographic 
effects  is  unknown.  The  maximum  surface  pressure  change  at  the 
ridge  is  -18  mb  h-^  during  the  switch-on  period.  A  quasi-steady 
flow  is  established  by  24  h,  with  surface  pressure  tendency  less 
than  0.3  mb  h-^  anywhere  in  the  model.  The  switch-on  of  the 
mountain  causes  a  surface  pressure  increase  of  ~1  mb  over  the  flat 
region.  A  similar  technique  to  Introduce  a  mountain  has  been 
utilized  in  Hayes  and  Williams  (1977)  and  Tibaldi  et  al  (1980). 

To  test  the  sensitivitiy  of  the  numerical  solution  to  the  moun¬ 
tain  switch-on,  we  have  conducted  several  supplementary  experiments 
with  a  uniform,  unsheared,  straight  basic  flow.  The  effects  of 
cumulus  convection  and  surface  friction  were  suppressed  for  these 
experiments.  First  we  varied  the  switch-on  period.  No  significant 
difference  was  found  in  the  quasi-steady  solutions  among  the  cases 
whose  switch-on  periods  vary  from  8-24  h.  We  also  varied  the  size 
and  height  of  the  topography.  It  was  found  the  details  of  the  flow 
were  sensitive  to  the  characteristics  of  the  mountain.  But  there 
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was  a  trend  In  cha  variations  and  tha  large  features  in  tha  flow 
ware  paraiatant.  Wa  thua  conclude  that  tha  idealised  topography 
and  tha  method  to  introduce  it  are  adequate. 


Tha  reaction  of  the  baaic  uniform,  unahaared  flow  to  tha 
idealized  topography  ia  of  couraa  vary  intareating  by  itaelf.  Since 
there  ia  no  vertical  wind  ahear  in  the  baaic  flow,  there  ia  little 
evidence  of  aountaln  wavea.  The  surface  flow  pattern  is  quite  different 
from  that  at  higher  levels.  The  surface  flow  is  highly  anticyclonic, 
with  much  of  the  mass  flowing  around  the  south  side  of  the  mountain  due 
to  its  blocking  effect.  There  is  a  wake  of  low  wind  speed  on  the  leeside 
until  the  northern  branch  and  the  stronger  southern  branch  converge  at 
250  km  west  of  the  ridge.  At  TOO  mb,  the  flow  is  still  anticyclonic  around 
the  mountain  with  a  leeside  trough.  There  is  little  branching  of  the 
flow  to  either  side  of  the  mountain.  Overall,  the  mountain-induced 
perturbation  at  TOO  mb  is  smaller  than  that  at  the  surface. 


4.  Orographic  Effects  on  the  Translating  Vortex 

Several  numerical  integrations  were  carried  out  to  isolate 
different  physical  processss.  Table  1  lists  all  the  axperiments. 
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Tab La  l 


List  of  All  Experiments 


Characteristics 

1 

Ho  topography 

2 

Topogrpahy 

l 

Topography,  Ho  Cumulus  Heating 

U 

Topography,  Ho  Surface  Friction 

Experiment  1  contains  no  topography.  Tho  tropical  cyclona  in  tha 
experiment  raachas  quasi-steady  stata  at  24  h  with  cantral  prassura 
of  978  nb  and  maximum  wind  of  32  •  a*  .  Its  path  la  approximately 
3°  to  tha  right  of  tha  naan  aaatarly  flow,  similar  to  the  result 
of  Chang  and  Madala  (i960}.  Tha  mean  translating  speed  is  that  of 
the  ambient  mean  flow.  In  this  section  ve  focus  on  the  comparison 
of  Sxps.  1  and  2. 

a.  Intensity  and  Path 

Tha  cantral  prassura  of  Exp.  2  is  never  lower  than  968  mb, 
which  is  10  mb  weaker  than  Exp.  1.  Tha  cantrsl  prassura  of  988  mb  is 
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reached  at  24  h  and  stays  with  ♦  2  tab  until  54  h  whan  the  Integration 
is  terminated.  Tha  maximum  surface  wind  la  greatly  affected  by 
the  mountain.  It  reaches  32  a  s”1  at  18  h,  decreasing  to  22  a  s 
Just  before  the  landfall  (Fig.  3).  But  as  the  center  of  the 
circulation  moves  past  the  topography,  the  storm  regains  It 
strength,  the  maximum  wind  increases  to  26  a  a"1.  The  decrease 
of  maximum  winds  as  typhoons  pass  Taiwan  Is  clearly  Illustrated  In 
Brand  and  Blelloch  (1974). 

The  mountain  has  great  Influence  on  the  translating  speed 
and  path  of  the  cyclone.  He  learned  that  the  tropical  cyclones  travel 
over  ocean  surface  with  a  mean  path  approximately  5°  to  the  right 
of  the  basic  currant  and  a  mean  speed  close  to  that  of  the  basic 
current  (Chang  and  Madala ,  1980;  also,  Exp.  1).  However,  In  Exp.  2, 
the  translation  speed  of  the  tropical  cyclone  Increases  as  It 
approaches  the  mountain  to  10  m  a**  (Fig.  4)  as  compared  to  the  speed 
of  the  basic  flow  of  5  a  s*1.  After  the  storm  center  passes  the 
ridge,  a  sudden  deceleration  takes  place.  By  54  h,  the  translation 
speed  decreases  to  -3  a  s'1. 

Because  of  the  positive  static  stability,  the  air  in  the  lower 
troposphere  tends  to  travel  around  the  mountains  Instead  of  over  It. 
When  the  storm  center  is  located  to  the  east  of  the  range,  a  stronger 
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•••early  flow  la  formed  to  th«  north  aid*  of  th«  mountain,  whereas 
the  westerly  wind  south  of  the  canter  is  turned  into  southerly. 

This  cyclonic  circulation  in  the  lower  troposphere  around  the 
mountain  not  only  speeds  up  the  storm  but  also  causes  the  storm  to 
move  north  to  the  east  of  the  mountain  (Fig.  6)  and  south  to  the  west 
of  the  mountain. 

The  changes  in  speed  and  the  curvature  in  the  path  of  the 
storm  are  in  agreement  with  observations  (cf.  Fig.  1). 

b.  Circulation 


The  wind  vector  and  vorticity  at  36,  42,  48  and  54  h  at  the 
700  mb  and  the  surface  levels  are  shown  in  Fig.  5  and  6.  At  the 
700  mb  level,  the  circular  wind  field  near  the  center  remains  nearly 
Intact  as  the  storm  passes  over  the  mountain  ridge.  The  easterly 
flow  north  of  the  island  is  much  stronger  than  the  flow  to  the  south, 
which  turns  from  westerly  at  36  h  to  southerly  at  48  b.  The  vorticity 
field  indicate  that  the  vorticity  is  in  general  cyclonic  in  the 
northern  part  of  the  island  mountain  because  of  its  proximity  to 
the  storm.  The  vorticity  is  in  general  cyclonic  everywhere  except 
on  the  windside  in  the  southern  part  of  the  mountain  where  it  is 
anti eye Ionic.  A  secondary  center  at  the  southern  tip  of  the  island 
exist 3  briefly  around  36  h  but  it  has  no  significant  consequences. 
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The  flow  at  th«  surface  level  (a  •  0.965)  is  more  complex. 

The  circular  wind  field  near  Che  storm  center  ie  deformed,  and  the 
vorticity  contours  stretch  along  the  ridge.  Due  to  the  orographic 
obstruction,  there  la  very  little  wind  component  normal  to  the  ridge, 
air  tends  to  flow  around  the  mountain.  This  results  in  a  sharp 
turning  of  airflow  on  the  wlndside  slope  and  a  wake  of  low  speed  on 
the  leeside  slope. 

At  42  h,  when  the  center  is  located  near  the  ridge,  the 
airflow  around  the  center  is  disorganized,  and  the  maximum  vorticity 
associated  with  the  storm  is  reduced.  Secondary  low  at  southern  tip 
of  the  mountain  exists  from  36-42  h.  Convergence  of  air  into  the 
main  center  is  evident  throughout  the  period.  At  54  h,  the  circular 
wind  field  near  the  center  is  reorganized,  and  the  vorticity  field 
had  regained  its  intensity. 

As  mentioned  before  the  path  of  the  storm  has  a  cyclonic 
curvature  induced  by  circulation  around  the  mountain.  The  storm's 
fast  translation  speed  between  36  and  48  h  is  also  due  to  its  northerly 
shift  into  the  stronger  easterly  to  the  north  of  the  mountain.  The 
circulation  of  the  surface  layer,  defined  as 


C  ■  &  u  dx  +  v  dy  , 


(4) 
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integrated  counterclockvisely  around  a  480  Van  x  720  tan  rectangular 
box  centered  around  point  D  in  Figure  2,  is  computed  (Fig.  7).  The 
surface  level  circulation  is  about  -4  x  10^  cm2s-1  at  24  h.  It 
rapidly  increased  with  time  to  +18  'u  20  x  10'*' °  cm2  s-1  at  40  hr» 
indicating  a  strong  cyclonic  flow  around  the  terrain. 


To  further  demonstrate  the  strong  easterly  around  the  north 
side,  the  rate  of  mass  flows  defined  as 

M(x)  -  -  i  f°'  f  P  (x,y)  u  (x,y,0)  dy  dCf  (5) 

g  1  y  8 

are  computed.  The  quantity  M  represents  the  rate  of  westward  mass 
flow  across  a  vertical  cross-section  at  certain  longitude  bounded 
by  a  *  0.5  and  1. 

Fig.  8  shows  the  values  of  M(x)  for  the  cross-section  from  the 
northerm  tip  of  the  island  to  the  northern  model  boundary  (solid 
line)  and  the  cross-section  from  the  southern  tip  of  the  island  to 
the  southern  model  boundary  (dashed  line).  Since  the  basic  flow  is 
uniform,  the  values  of  the  two  M(x)  are  equal  at  -v  2.8  x  1013  g 
along  the  lateral  inflow  boundary.  During  the  period  24  h  *v  48  h, 
the  northern  M(x)  increases  to  5  ^  6  x  10  g,  while  the  southern 
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M(x)  decreases  to  1  x  10  g  or  even  negative  values  Indicating  a 
reversed  (eastward)  mass  flow  there. 

It  Is  evident  from  these  results  that  the  mountain  has 
significant  Influence  on  the  propagating  tropical  cyclones.  The 
terrain  induces  a  low-level  cyclonic  circulation  around  the  mountain, 
and  diverts  most  of  the  easterly  mass  flow  to  the  north  of  the 
terrain.  Thus  the  storm  center  is  forced  to  move  cyclonically 
about  the  mountain  into  the  strong  easterly  flow  there  and  results 
in  its  fast  movement  across  the  mountain.  The  mountain  weakens  the 
circulation  by  practically  blocking  the  airflow  at  the  low  levels. 

c.  Vertical  Velocity  and  Precipitation 

The  topography  has  strong  effects  on  the  vertical  velocity 
field,  the  low-level  upward  motion  mainly  occurs  on  the  wlndside 
slopes.  As  evidence  by  the  u>(-dp/dt)  field  at  a  *  0.97  (Fig.  9), 
the  peak  upward  motion  at  36  h  occurs  in  the  region  to  the  north  of 
the  storm  center  and  over  the  northern  part  of  the  island,  where  the 
wind  direction  is  in  land  (see  Fig.  6).  At  42  h,  the  region  with  the 
upward  motion  is  widespread,  but  mainly  restricted  to  the  northern- 
part  of  the  mountain  range  where  wind  speed  is  high.  At  U8  h,  the  strong 
upward  motion  is  shifted  to  the  windslopes  in  the  south  of  the.  mountain. 
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The  upward  motion  at  PBL  top  Is  well  correlated  with  precipi¬ 
tation,  which  is  responsible  for  most  of  the  typhoon  damage.  The  real 
precipitation  distribution  of  course  is  more  complicated  due  to  the 
complex  terrain  and  the  large  scale  distribution  of  specific  humidity.  The 
distributions  shown  in  Fig.  9  construct  only  a  general  picture,  because 
the  initial  specific  humidity  field  in  the  model  is  assumed  uniform. 


5.  The  Effect  of  Latent  Heat  Released  in  Connective  Clouds 


To  determine  the  effect  of  the  release  of  latent  heating  by 
the  cumulus  convection,  we  Integrated  the  numerical  model  with  the 
diabetic  heating  and  the  surface  friction  terms  suppressed  (Exp.  3). 
Without  there  two  effects,  we  effectively  eliminate  the  possibility 
for  the  initial  vortex  to  intensify.  Indeed,  the  vortex  experiences 
a  gradual  weakening  throughout  the  integration  due  to  internal  friction. 

The  movement  of  the  low-level  vortex  center  and  the  airflow 
around  the  mountain  in  Exp.  3  are  notably  different  from  those  in 
Exp.  2.  The  original  westbound  surface  vortex  center  (Fig.  10) 
decelerates  between  36  and  48  h  as  it  approaches  the  mountain  and 
shifts  northward.  It  stops  at  the  location  denoted  by  60  h  east  of 
the  mountain.  It  dissipates  completely  by  72  h.  As  the  original  vortex 
is  dissipating,  a  surface  secondary  center  appears  at  36  h  at  the  southern 
tip  of  the  mountain  and  moves  slowly  northward. 
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At  60  h,  another  surface  secondary  center  occurs  300  km  to  the  west. 

In  the  mean  time,  the  700  mb  center  passes  the  ridge  at  the  speed  of 
the  basic  current,  taking  a  slightly  antlcyclonic  trajectory.  The 
secondary  center  develops  and  becomes  the  major  surface  center  as  the 
upper  level  center  aligns  vertically  with  it  over  the  ’open  sea  after  72  h. 

There  Is  a  strong  easterly  airflow  with  antlcyclonic  vortlclty 
to  the  south  of  the  mountains  In  the  surface  flow  (Fig.  11). 

This  antlcyclonic  low  with  minimum  vortlclty  to  the  southeast  of  the 
mountain,  extends  northward,  deforms  the  vortex,  and  shifts  the 
vortex  center  northward.  The  wake  to  the  west  of  the  ridge, 
characterized  by  low  wind  speeds,  contains  return  flows  and  secondary 
center.  The  upper  airflows  (not  shown)  are  similar  to  those  In 
Exp.  2  (Fig.  5)  except  for  the  antlcyclonic  curvature  in  the  path  of 
the  vortex. 

The  low-level  circulation  around  't^e  terrain  (computed  by 
Eq.  4,  plotted  in  Fig.  7)  is  predominantly  antlcyclonic  at  -  8  *  -  10  x 
10  ^  cm^  s  \  except  for  a  brief  period  -  60  h.  Even  then,  the 
maximum  value  Is  no  more  than  2  x  10*^  cm^  9  \  The  values  of  M(x) 
(defined  by  Eq.  3,  plotted  in  Fig  8)  show  much  less  north-to-south 
polarization  compared  with  Exp.  2  due  to  the  presence  of  the  stronger 
easterly  surface  flow  to  the  south  and  the  absence  of  the  strong 
easterly  flow  to  the  north. 
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The  momentum  field  in  Exp.  3  i3  decoupled  vertically,  the  wind  shear 
between  the  surface  level  and  upper  level  is  unrealistically  large. 

Exp.  3  of  course  cannot  represent  any  real  situation  because  there  is 
always  vertical  coupling  due  to  convection  and  turbulence.  However, 
according  to  Wang's  (1980)  classification,  there  are  typhoons  that  passed 
Taiwan's  CMR  (Type  I,  several  examples  are  shown  in  Fig.  1)  and  there 
are  those  that  were  blocked  by  CMR  and  reformed  from  secondary  center 
(Type  II).  Figure  12  shows  some  type  II  cases.  Table  2  shows  that  the 
intensities  of  types  I  and  II  before  the  landfalls  are  quite  different. 

The  mean  central  pressure  of  type  I  is  30  mb  lower  than  type  II,  the 
mean  maximum  wind  of  type  I  is  M.8  m  s”1  stronger  than  type  II.  The 
observed  typhoons  are  more  intense  than  our  simulation  mainly  because 
of  the  horizontal  resolution  in  our  model.  We  have  performed  further 
tests  on  the  resolution  problem  and  the  results  are  presented  in  the  Appendix. 

It  is  conceivable  that  there  are  stronger  cumulus  convection  and 
better  vertical  structural  coherence  in  strong  typhoons  than  weak 
typhoons.  Our  experiments  suggest  that  intensity  is  a  deciding  factor 
in  tropical  cyclone's  behavior  in  passing  high  mountain  range.  Strong 
cyclones  pass  with  modifications  and  weak  cyclones  are  disturbed 
to  a  large  extend  with  their  low-level  circulations  blocked. 
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TABLE  2 


Intensities  of  Two  Types  of  Typhoons  Passing  Taiwan 


Type 

Number  of 

Cases 

Mean  Central  Pressure 

(Standard  Deviation) 

Mean  Maximum  Wind 

(Standard  Deviation) 

-1 

I 

22 

949  mb 

52.4  m  s 

(23.8) 

(9.0) 

II 

12 

979  mb 

36.8  m  s 

(12.7) 

(12.6) 

We  thus  conclude  that  the  latent  heating  by  cumulus 
convection  plays  an  Important  role  in  the  tropical  cyclone's 


reaction  to  the  topography.  The  latent  heating  keeps  the 
circular  wind  field  relatively  intact  in  spite  of  the  orographic 
obstruction.  It  produces  a  stronger  easterly  flow  to  the  right 
side  (relative  to  the  storm  track)  of  the  mountain  with  cyclonic 
circulation.  Such  airflow  field  produces  a  fast  translation  speed 
and  a  cyclonic  movement  of  the  storm  center  as  shown  by  Fig.  1. 

The  lack  of  latent  heating  results  in  the  weakening  of  the  storm. 
When  encountering  the  mountain  range,  weak  tropical  cyclones  tend 
to  form  low-level  anticyclonic  circulation  around  the  mountain  and 
develop  leeside  secondary  center. 

6.  Effect  of  Surface  Friction 


Exp.  4  is  conducted  with  the  cumulus  parameterization 
suppressed.  A  comparison  of  Exps.  3  and  4  shall  indicate  the 
effect  of  the  surface  friction.  The  behavior  the  vortex  in 
Exp.  4  is  very  similar  to  that  in  Exp.  3  in  the  vortex's 
translation  speed,  path  and  flow  pattern  at  high  levels.  The 
surface  flow  in  Exp.  4  is  weaker  than  that  in  Exp.  3  as 
expected.  The  original  vortex  is  hardly  Identifiable  at  48  h, 
the  secondary  center  is  also  less  well  defined.  The  weak 
circulation  is  illustrated  by  Fig.  7. 
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We  should  note  that  Exp.  4  again  does  not  represent  any 
realistic  situation.  The  surface  friction,  which  dissipates 
kinetic  energy,  produces  energy-generating  cross-isobaric  flow. 
The  frictionally-induced  convergency  further  enhances  the 
cumulus  convection  in  conditionally  unstable  atmosphere. 

However,  Exp.  4  points  out  that  the  surface  friction,  over 
either  the  mountain  or  the  ocean,  does  not  contribute  directly  to 
the  particular  behavior  of  a  tropical  cyclone  passing  over  a 
mountain. 


Vorticity  Budget 


It  is  useful  for  field  forecasters  to  estimate  from  a 
simplified  vorticity  budget  the  circulation  (or  mean  vorticity) 
around  the  island  mountain  when  a  typhoon  approaches.  Neglecting 
the  dissipation  for  the  time  being,  one  can  approximate  the 
absolute  vorticity  tendency  following  an  air  volume  by 


~  (C  +  f)  -  -  (C  +  f)  V  •  V  -  -  (;  +  f)  (V  •  V)  -  (5  +  f>  CV  *  V)  (6) 
dt  uo,  U  U  °  u-ud 


where  (V  *  V)  is  the  orographically-induced  divergence  and 
'u  a,  o 

*  V)  is  the  divergence  generated  by  diabatic  processes  in 
the  storm.  To  maintain  the  obsolute  vorticity  in  the  air  volume, 
one  needs 


“(V  *  V)  >  (V  *  V> 

%  \  d  “  %  %  o 


(7) 
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Su... 


Suppose  air  flows  completely  over  the  mountain  at  raid-troposphere 
the  orographically-induced  divergency  on  the  vindside  slope  can 
be  expressed  for  first-order  approximation  as 

V  9<j) 

(V  .  y)  _Q_  ( - 2-)  (8) 

gH  9x 


where  Vn  is  a  typical  wind  speed  normal  to  the  mountail  range  near 

1  3$ 

a  tropical  cyclone,  H  is  the  scale  height.  Expression  -  ( — -&-) 

g  a* 

is  the  slope  of  the  mountain,  and  is  equal  to  0.0167  for 

the  idealized  topography.  Let  Vn  ■  25  m  s  ^  and  H  *  8  km,  we  get 
-  (V  •  V)  £  -5  x  10  S1,  which  is  equivalent  to  a  cumulus 

heating  rate  of  ~  110  °K  day  or  precipitation  rate  of 
,,  35  cm  day  This  means  that  the  divergence  induced  by  our 

idealized  topography  is  equivalent  in  magnitude  to  the  convergence 
in  a  fairly  strong  tropical  system.  In  a  quasi-steady  tropical 
cyclone,  the  vorticity  generated  by  -  (V  *  V),  is  nearly  balanced 

v  'u  a 

by  internal  dissipation  and  surface  friction.  By  obtaining  (V  *  v) 

'b  'bo 

from  (8),  one  can  estimate  the  rate  of  exponential  decay  in  time  of 
storm's  vorticity  by  integrating  a  simplified  vorticity  equation. 


In  reality,  the  vorticity  budget  is  more  complex  than  the 
disscusion  above.  Having  a  three-dimensional  model,  we  can 
diagnose  the  contribution  to  vorticity  generation  of  each  term  in 


the  vorticity  equation 


The  vortieity  equation  in  pressure-coordinates  can  be  written 
as  (Holton,  1979) 


3£  .  «  .  7  (r  +  f)  -  w  3?  -  (C  ♦  f)  V  •  v  +  k  .  (W  x  +  S  (9) 

3t  a*  '  3p  ^  ^  3p 

(I)  (II)  (III)  (IV)  (V) 

where  symbols  have  their  conventional  meterorologlcal  meanings. 

The  terms  on  the  right  represents,  respectively,  the  effects  of 

horizontal  advection,  vertical  advection,  divergence,  twisting, 

and  sink  or  source  due  to  diabatic  processes.  The  last  term 

is  mainly  due  to  friction  and  diffusion  in  our  model,  it  should 

be  proportional  to  the  value  of  the  vortieity  due  to  the  physical 

parameterization  in  our  model.  The  characteristic  scales  pertinent  for 

our  present  study  are  U.-.5  m  s  \  W  10  mb  h  \  L  ^  200  km,  Ap  t 

100  mb  and  f  *  5  x  10  ^  s  \  A  straightforward  scale  analysis 

-9 

gives  orders  of  magnitude  for  the  terms:  term  I  -  0  (10  ), 

term  II  -  0  (10_9),  term  III  ~  0  (10-8),  and  term  IV  -  0  (10-9). 

The  magnitude  of  term  V  is  difficult  to  estimate  for  the 
problem  under  investigation,  it  should  be  equal  to  or  less  than  the 
magnitude  of  term  I  or  III.  Therefore,  unlike  the  large  scale  flow, 
the  effects  of  twisting  and  divergence  are  far  from  negligible  in 
the  vortieity  budget  for  this  mesoscale  problem. 
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Because  Che  cyclonic  circulation  around  Che  mountain 
grows  rapidly  between  24  to  48  h  (Tig.  7)  and  the  surface  pressure 
at  the  peak  is  about  800  ob,  we  selected  the  700  mb  flow  at 
36  h  (Fig.  5)  for  vorticity  budget  study.  In  order  to  make 
use  of  Eq.  (9),  model  data  are  linearly  interpolated  to  isobarlc 
surfaces.  Vertical  differences  are  taken  over  a  ISO  mb  thickness. 

It  turns  out  that  the  contribution  of  vertical  advectlon  is  small. 
The  distributions  of  terms  I,  III  and  IV  are  shown  in  Fig.  13. 

We  note  that  the  36  h  vorticity  field  (Fig.  5a)  consists  of 
strong  cyclonic  vorticity  near  the  center  and  to  the  leeside 
slopes,  and  anticyclonic  vorticity  on  the  windside  slopes.  The 
contribution  to  local  vorticity  change  due  to  horizontal  advectlon 
(Fig.  13a)  is  approximately  1/4  wavelength  out  of  phase  downwind 
from  the  vorticity  field.  The  strong  positive  tendency  between  the 
ridge  and  the  storm  center  in  conjunction  with  the  positive  tendency 
downwind  from  the  leeside  trough,  contributed  to  the  formation  of 
the  discussed  cyclonic  circulation  around  the  mountain  (Fig.  7) 
at  this  time. 

There  is  negative  tendency  in  the  distribution  of  term  III 
(Fig.  13b)  over  the  windside  slopes  and  positive  tendency  over  the 
leeside  slopes  due  to  vortex  stretching.  The  positive  tendency 
northwest  to  the  storm  center  is  apparently  associated  with  the 
tropical  cyclone  itself.  The  high  wind  speed  and  cumulus  con¬ 
vections  in  the  front-and-right  quadrant  produces  a  strong 
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convergence  there,  end  hence  the  positive  tendency.  This  positive 
tendency  maintains  s  strong  vorticity  neer  the  store  center,  which 
in  turn  produces  the  distribution  of  term  I  due  to  horizontal 
advection.  In  Exp.  3,  where  cumulus  convection  and  surface 
frictin  is  suppressed,  this  legion  of  positive  tendency  is 
conspicuously  absent.  Note  also  that  the  center  of  the  positive 
tendency  region  is  located  north  of  the  storm  -.iter,  this  is  a 
symptom  of  the  northward  movement  of  the  center.  There  is  another 
maximum  positive  tendency  to  the  south  of  the  mountain  associated 
with  the  secondary  center  (see  Fig.  6a). 

The  twisting  term  (term  IV)  is  confined  to  the  mountainous 
region.  It  indicates  that  only  in  this  region  a  horizontal 
differential  of  vertical  motion  is  coupled  with  a  vertical  wind 
shear  and  that  the  flow  near  the  mountain  is  highly  baroclinic. 

Again  we  find  in  the  vorticity  budget  that  the  dlabatic 
effects  Play  an  important  role  in  producing  the  orographic  effect. 
The  diabetic  processes  generate  convergence  near  the  center  of  the 
tropical  cyclone  to  maintain  its  vorticity  field,  which  by 
horizontal  advection,  causes  a  cyclonic  circulation  around  the 
island  mountain. 
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8.  Summary  and  Discussion 

The  orography  effects  on'  propagating  tropical  cyclone 
exerted  by  an  Idealized  mesoscale ,  Island  mountain  range  was 
Investigated  by  using  a  multi-layer^-PB-  model  with  parameterized  physics. 
The  numerical  simulation  successfully  reporduces  the  following 
observed  behavior  of  the  westward  propagation  tropical  cyclones: 

A' 

(1)  The  Intensity  of  the  tropical  cyclone  is  reduced;  (2)  the 
translating  speed  increases  before  passing  the  mountain  ridge; 

(3)  the  storm  center  moves  cyelonically  in  passing  over  the  mountain; 

(4)  the  distribution  of  vertical  motion  is  orographically 
modulated,  upward  motion  occurs  on  the  windside  slopes  where 
the  vorticity  is  positive. 

The  results  indicated  that  the  approaching  tropical 
cyclone  induces  a  cyclonic  circulation  around  the  topography. 

The  mountain  range  constitutes  a  physical  obstruction  to  the 
low-level  flow  and  generates  a  strong  easterly  mass  flow  to  the 
north  of  the  mountain  range,  r- 

Further  experiments  showed  that  the  strong  latent 
heating  by  cumulus  convection  is  responsible  for  the  generation 
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of  the  above  mentioned  cyclonic  circulation.  In  the  experiment 
without  cumulus  heating,  the  circulation  around  the  mountain  is 
mostly  anticydonlc,  and  the  low-level  vortex  center  is  blocked 
by  the  mountain  range.  Secondary  vortex  centers  form  in  the  lee 
trough  and  develop  when  in  phase  with  the  upper  level  vortex 
center.  These  experiments  helped  explain  the  observed  different 
behavior  of  weak  typhoons  from  that  of  stronger  ones. 

The  important  role  of  the  dlabatlc  effects  was  confirmed 
by  a  vorticity  budget  study  on  the  flow  before  the  landfall. 

It  is  the  convergence  generated  by  the  cumulus  heating  and  surface 
friction  that  produces  a  strong  vorticity  field.  The  positive 
vorticity  advection  in  front  of  the  storm  compensates  the  negative 
vorticity  tendency  due  to  vortex  compression  on  the  windside 
slope,  and,  in  conjunction  with  the  vortex  stretching  on  the 
leeslde,  results  in  a  positive  mean  vorticity  over  the  mountain 
range. 


Our  numerical  model  can  of  course  be  improved.  The 
major  shortcomings  are  related  to  the  PBL  formulation  over  the 
mountain  and  the  horizontal  resolution.  We  may  have  stretched 
the  model’s  PBL  parameterization  in  applying  it  over  rugged 
mountain  terrains  under  windy  conditions.  A  multilevel  PBL 


28 


formulation  is  perhaps  needed  in  spite  of  the  cost.  The  60  km 
resolution  in  our  model  is  only  marginally  acceptable  for  the 
storm  (see  Appendix).  A  nested  grid  system  can  help  circumvent 
the  problem  at  not  much  greater  expense. 

Many  of  our  conclusions  cannot  be  substantiated  due  to  a 
lack  of  indepth  analyses  of  observation  data.  Case  studies  as 
well  as  composite  studies  on  existing  data  are  needed  to  further 
advance  our  understanding  of  the  tropical  cyclone's  behavior  in 
crossing  mountain  ranges. 
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Fig.  1  —  The  Terrain  of  Taiwan  and  paths  of  some  typhoons. 
The  contours,  plotted  every  1000  m,  are  simplified. 


SO  KM 


Fig.  2  —  The  idealized  topography:  (a)  a  plan  view,  shaded  regions  are  elevated; 
the  dark  line  in  the  center  represents  the  ridge;  dots  are  model  grid  points; 
(b)  a  northward  view  of  the  AB  cross-section;  (c)  an  eastward  view  of  the  CD 
cross-section. 


12  18  24  30  36  42  48  54 

t(h) 


Fig.  3  —  The  temporal  variations  of  the  maximum  surface  wind  of 
Exp.  2.  The  shaded  area  denotes  the  time  the  strom  center  is 
over  the  land. 


Fig.  4  —  Same  as  Fig.  3  except  for  the  6-h  average  translating  speed. 
Note  the  speed  of  basic  flow  is  ~  5  ms  _1. 
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Fig.  5  —  The  700-mb  wind  vectors,  relative  vorticity  field,  present  and  past  center 
positions  of  Exp.  2  at  (a)  36  h,  (b)  42  h,  (c)  48  h,  and  (d)  64  h.  The  values  of  the 
positive  vorticities  (solid)  are  1,  2,  4, ...  X  10~4  s"1.  The  value  of  the  negative  vor- 
ticities  (dash)  are  -1,  -2,  -3, ...  X  10-4  s'1. 
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Fig.  6  —  Same  as  Fig.  5  except  for  the  surface  level  ( a  *  0.965) 
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Fig.  7  —  The  surface  level  circulation  along  a  480  X  640  km  rectangular  box 
centered  around  the  ridge  for  Exps.  2,  3,  and  4 


Fig.  8  —  The  man  flow  as  defined  by  Eq.  (5)  of  the  northern  (solid)  cross-section 
and  southern  (dash)  cross-section  for  Exps.  3  and  3.  Positive  values  denote  cy¬ 
clonic  circulation;  negative  anticydonic. 
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Fig.  9  —  The  vertical  velocity  at  a  ■  0.93  at  (a)  36  h,  (b)  42  h,  and  (c)  48  h.  Contour 
intervals  is  25  mb  h-1 .  Areas  of  upward  motion  stronger  than  -25  mb  h_1  are  shaded. 
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Fig.  10  —  The  movement  of  the  surface  ($)  and  700-mb  (o)  centers  for  Exp.  3.  The 
movements  of  secondary  centers  ($)  are  traced  with  a  dashed  line. 
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Fig.  11  —  The  surface  wind  vector  and  vorticity  field  at  48  h  for  Exp.  3  showing  the 
deformation  of  the  vortex  and  the  strong  anticycionic  flow  around  south  of  the 
mountain 
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Fig.  12  —  Paths  of  some  weak  typhoons  that  were  prevented  from  passing  the 
Central  Mountain  Range,  and  were  replaced  by  the  leeside  secondary  centers 
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Fig.  13  —  The  spatial  distributions  of  (a)  horizontal  advection,  (b)  divergence,  and  (c)  twisting 
effects  in  the  vorticity  budget  at  700-mb  at  36  h  for  Exp.  2.  Contour  interval  is  5  X  10~9  s’2. 
Solid  lines  for  positive  values,  dashed  line  for  negative  values. 


APPENDIX 


It  is  generally  accepted  that  the  intensity  of  simulated 
tropical  cyclones  depends  on  the  horizontal  resolutions  of  the 
numerical  models.  High  horizontal  resolutions  near  the  center 
are  necessary  to  resolve  the  eye  structure  and  strong  pressure 
gradient  there.  The  high  angular  momentum  in  the  inflow  can  also  be 
realized  at  small  radii  in  models  with  high  horizontal  resolutions. 
Rosenthal  (1970)  found  differences  in  storm  structure,  kinetic 
energy  production,  and  rainfall  rate  even  between  a  model  with 
20-km  resolution  and  another  with  10-km  resolution. 

In  order  to  show  that  the  difference  in  the  intensity 
between  our  model  tropical  cyclone  and  that  observed  (Table  2) 
is  related  to  the  60-km  resolution  in  our  model,  a  set  of  supple¬ 
mentary  numerical  experiments  were  carried  out.  The  model 
used  is  an  axisymmetric  analog  to  our  three-dimensional  model. 

The  initial  condition,  environmental  stability,  sea-surface 
temperature,  model  physics  are  kept  the  same  as  the  three- 
dimensional  model.  Three  integrations  were  performed  with  a  15-km 
(Exp.  A15),  a  30-km  (Exp.  A30),  and  a  60-ka  (Exp.  A6o)  resolution. 
Quasi-steady  states  were  reached  in  both  experiments  after  30  h. 

A  comparison  of  important  storm  parameters  at  36  h  is  listed  in 
Table  Al. 
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TABLE  A1 


r 


Parameters 

A15 

A30 

A60 

Min.  Pressure  (mb) 

9^1 

9U5 

973 

-1 

66 

Max.  Wind  Speed  (ms  ) 

57 

34 

Radius  of  Max.  Wind  (km) 

1*5 

60 

120 

Max.  Rainfall  Rate  (cm  day”^-) 

185 

167 

59 

-1 

-65 

-U6 

Max.  Upward  Velocity  (cb  h  ) 

-19 

Average  Rainfall  Within 

r  ■  120  km  Ccm  day  S 

60 

53 

16 

Conversion  Rate  from  P.E.  to 

K.E.  within  r  -  300  km  (1012  W) 

u.o 

5.5 

3.3 

There  is  clearly  considerable  differences  among  the  three 
integrations.  But  the  difference  between  A15  and  A30  is  smaller 
than  that  between  A30  and  a60.  The  stormy  intensities  of  A15  and 
A30  are  similar  to  those  observed  in  typhoons  (Wang,  1980),  whereas 
the  intensity  of  A60  is  similar  to  our  Exp.  2,  The  rate  of  P.  E. 
conversion  to  K.  E.  in  A15  within  300-km  radius  is  smaller,  because 
the  storm  in  A15  is  more  concentrated. 

These  values  suggest  that  a  three-dimensional  model  with  15-km 
or  30-km  resolution  can  produce  a  storm  intensity  comparable  to  those 
observed  (Table  2)  and  that  the  coarse  resolution  in  our  three-dimensional 
model  is  responsible  for  the  less  intense  tropical  cyclone  in  Exp.  2. 

It  is  perhaps  advisable  to  use  a  30-km  resolution  in  our  three- 
dimensional  model,  but  economic  considerations  dictated  the  60-km 
resolution. 
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